Near Equipartition Magnetic Fields Measured in the Cool Cores of Galaxy Clusters 
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Tangential discontinuities, seen as X-ray edges known as cold fronts (CFs), are ubiquitous in 
cool-core galaxy clusters. We analyze all 18 deprojected CF thermal profiles found in the literature, 
including three new CFs we identify (in clusters A2204 and 2A0335). We discover small but sig- 
nificant thermal pressure drops below all non-merger CFs, and argue that they arise from strong, 
near-equipartition (10-20%) magnetic fields below and parallel to the discontinuity. Such mag- 
netization can stabilize the CF against Kelvin-Helmholtz instabilities, and explain the connection 
between CFs and radio minihalos. 



In the past decade, high resolution X-ray observations 
have revealed an abundance of density and temperature 
discontinuities known as cold fronts (CFs). They are 
broadly classified as cool-core CFs and merger CFs (for a 
review, see [H); a third, putative class of shock-collision 
induced CFs [2J, has not yet been observed. Here we 
focus exclusively on the first, cool-core, type of CFs. 

Core CFs are observed in most of the otherwise re- 
laxed, cool-core clusters (CCs) 0, in distances up to 
~ 400 kpc from the center. They are usually nearly con- 
centric or spiral, and multiple CFs are often observed in 
the same cluster. The temperature contrast across such 
a CF is T /Ti ~ 2 where inside/outside subscripts i/o 
refer to regions closer to/farther from the cluster center, 
or equivalently below/above the CF. The plasma beneath 
the CF is typically denser, colder, much lower in entropy, 
and higher in metallicity, than the plasma above it. 

Such CFs are thought to be a quasi-spiral tangential 
discontinuity surfaces seen in projection [5|,l6|. They may 
reflect large-scale "sloshing" oscillations of the intraclus- 
ter medium (ICM), driven by mergers hi, possibly in- 
volving only a dark matter subhalo @, Q , or feedback 
from the central active galactic nucleus (AGN) in the 
form of weak shocks or acoustic waves displacing cold 
central plasma 0, They were also proposed to be 
the signature of long-lived spiral bulk flows underlying 
cool cores, sustained for example by AGN bubbles |6J. 

Deprojected thermal profiles across core CFs reveal the 
presence of fast, nearly sonic flows beneath, ami tangen- 
tial shear layers extending below, the CFs 
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shear can produce the near-equipartition magnetic fields 
needed to stabilize the CF against Kelvin-Helmholtz in- 
stabilities (KHI), rj B = P B /Pth = B 2 /8irP th = 10-20% 
[H) . Here, B is the magnetic amplitude, and equiparti- 
tion is defined as equal magnetic and thermal pressures, 
Pb = Pth- Such shear layers and magnetization levels 
were reproduced in CC sloshing simulations (l2j |. 

There is circumstantial evidence for strong magnetic 
fields along and below CFs. The thinness of the emission 
discontinuity, which is much narrower than the Coulomb 
mean free path, requires a magnetic suppression of trans- 
port across the CF 13]. The CF-radio minihalo connec- 



tion (Q; see also [H, EH), m particular when combined 
with the correlation between the radio and X-ray bright- 
ness of mini-halos, suggests strong magnetization below 
CFs, significantly exceedin g th e cosmic microwave back- 
ground equivalent ~ 3 /zG |17| . 

Previous measurements of core magnetic fields [Til - [2lj | 
produced conflicting results, and are not particularly sen- 
sitive to the CF region. In this letter we report the dis- 
covery and direct measurement of a thermal pressure Pth 
discontinuity across core CFs. We interpret this as a dis- 
continuous jump in nonthermal pressure P nt below the 
CF, and argue that it must be predominantly magnetic. 
We assume a Hubble constant of H = 70 km s _1 Mpc -1 , 
and a 0.23 He mass fraction. Error bars represent la con- 
fidence levels. 

Core CFs do show a thermal pressure discontinuity. — 
Unlike merger CFs, core CFs were so far thought to in- 
volve no change in thermal pressure (i.e. appear ther- 
mally isobaric; [l| and references therein). To critically 
examine this, we analyze the thermal profiles near core 
CFs, and derive the pressure on both sides of the CF. 

We extract all the thermal profiles across core CFs in 
the literature, and select only profiles that have been de- 
projected along the line of sight. There are 15 such CFs 
from 10 CCs, in which P t h can be estimated on both sides 
of the CF. We compute P t h as the product of the de- 
projected particle number density and temperature, as- 
suming a single-phase ideal gas. (Deviations from this 
assumption are discussed towards the end.) Some Pth( r ) 
profiles are illustrated in Figure [3] 

In FigureQ]we present the ratio £ = Pi/P between P t h 
just below the CF (Pi) and P t h just above the CF (P ), 
for all CFs, as a function of the CF radius r c f normalized 
by i?5oo ■ The values of -R500 , the radius enclosing a mean 
density 500 times the critical density of the Universe, are 
taken from Ref. [22J. Literature references for the CF 
data are provided in the caption. 

We first compute Pi and P assuming a radial power- 
law behavior on each side of the CF, estimated from the 
two nearest data points. In the three cases where only 
one temperature data point exists below the CF (the two 
CFs in RXJ1347.5 and the inner CF in A1644, marked 
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FIG. 1: The deprojected thermal pressure ratio £ = Pi/P , 
plotted against the normalized radius r/Rsoo, for all core CFs. 
The solid line corresponds to thermally isobaric, £ = 1 tran- 
sitions. Horizontal lines show the averages of all CFs (dotted 
line), of our final CF sample (dashed, Eq. Q), and of the fic- 
titious sample (dot-dashed); see text. An isothermal profile is 
assumed for CFs marked by empty symbols. Data: Ref. [25 1 
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by empty symbols in Figure [1} , we assume an isothermal 
temperature profile. We also identify three new CFs, seen 
coincident ally in the emission, temperature, and some- 
times also metallicity: two in the Southwest sector of 
2A0335 [23j, at radii ~ 28 and ~ 62 kpc, and one in the 
West sector of A2204 [H, at r cf ~ 127 kpc. These CFs 
are shown as dashed rectangles in Figure 03 

Confidence levels are estimated by error propagating 
the quoted uncertainties in the temperature and density 
data, assumed for simplicity to be statistical, normally 
distributed [is} , and uncorrelated. The resulting error is 
dominated by the uncertainty in the temperature bin just 
above the CF. Our error estimate is conservative, yielding 
high uncertainties in the inferred pressure jumps. A more 
realistic estimate, using the covariance matrix between 
points above and below the CF, should yield the same 
nominal points with a reduced uncertainty. 

All CFs, except two unusual cases discussed below, 
show £ < 1, suggesting a P t h deficit below the CF. The 
statistical significance per CF is modest (between 0.3<r 
and 3.6a, except Virgo with 7.3ct), but the similar trend 
among CFs renders this result robust. The mean result 
for all CFs is 1 all = 0.78 ± 0.03. The two unusual CFs 
with f > 1 are the outer CFs in A1644 and RXJ1347.5, 
both of which are known to be undergoing a merger. This 
suggests that these CFs were disrupted by the merger; 



excluding them from the analysis leaves a non-disrupted 
CF average J = 0.74 ± 0.03. The CF sector in Virgo is 
narrow and shows evidence of density substructure, as 
seen in Figures [2] and G2 so it is henceforth excluded as 
well. The mean result for the remaining 15 CFs is 



£ = 0.82 ±0.04. 



(1) 



We conclude that £ ~ 0.8, with a large dispersion among 
CFs, is significantly smaller than unity. 

In order to show that the result is not an artifact of the 
deprojection method, we apply the same analysis method 
to the projected radial profiles. This is possible in the 
clusters A496 and A133, for which the projected temper- 

n 



The resulting 



ature profiles are also published 
pressure profile are smoother due to the projection, but 
the pressure discontinuity, albeit reduced, is still seen. 

To examine the sensitivity of the results to our two- 
point power-law fit, we examine a different fit procedure. 
Instead of using only the two nearest data points on each 
side of the CF, we include all points of radius within a 
factor 2 from the CF. We avoid any additional CF, so 
a single power-law may still provide a good approxima- 
tion. All results change within their computed error bars, 
with a slight tendency towards £ values closer to, but still 
smaller than, unity. The average in this case increases, 
due to the smoothing effect of the fit, to £ = 0.85 ± 0.04. 

As a general test against possible systematic errors (in- 
cluding the above deprojection and fit errors), we com- 
pute £ with the same data and in the same method used 
to produce Figure [TJ but for radii where no CF is present. 
This is possible in five CC sectors, in which we place a 
fictitious CF between every two consecutive radial bins 
away from the true CFs. The resulting £ values, shown 
in Figure [2j lie above the CF mean of Eq. (CQ) for all but 
two fictitious CFs (both in the narrow sector in Virgo, ex- 
cluded from the above analysis due to suspected density 
substructure). The mean value here is £ = 1.05 ± 0.03; 
without Virgo it becomes £ = 1.03 ± 0.06. We simi- 
larly analyze the d eproje cted, azimuthally-averaged ra- 
dial profiles (from 32h34| ) . which are sensitive to some 



systematic errors [50|. This gives £ = 1.04 ±0.01. 

We conclude that there are small but significant drops 
in P t h below CFs, given on average by Eq. (CQ). None of 
the systematic effects considered can mimic this result. 

The missing pressure is nonthermal. — A drop in the 
total pressure below the CF would imply some radial mo- 
tion, perpendicular to the CF surface. As the flow must 
parallel the tangential discontinuity, such a drop would 
correspond to inward acceleration of the CF, or outward 
ram pressure of the gas below the CF [l[. However, a 20% 
pressure jump would correspond [l]| to Mach M. ~ 0.5 
radial motion: unrealistic in an otherwise relaxed core. 
Moreover, core CFs are thought to be part of extended 
three-dimensional surfaces spanning much of the core, so 
the implied coherent radial motion is unlikely. 
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FIG. 2: Fictitious CFs. Horizontal £ > 1 lines show the aver- 
ages £ of these CFs, with (dotted) and without (dot-dashed) 
Virgo. The other lines are defined in Figure [T] 



Hence, the total pressure must be nearly continuous 
across the CF, so the missing pressure must arise from 
some nonthermal component. The pressure drop we find 
imposes a lower limit on the nonthermal pressure P n t 
below the CF, as our measurement is not sensitive to any 
additional, smoothly varying nonthermal component. 

There is prior evidence for nonthermal pressure in clus- 
ters. Comparing optical and X-ray data gives r\nt = 
Pnt/Pth < 10% m the cores of Virgo and Fornax [35j. 
Comparing X-ray and weak lensing data yields r\ nt ~ 30% 
in MS2137 [36j|. However, these estimates assume hydro- 
static equilibrium, and are not specific to CFs. 

The distribution of nonthermal pressure. — To further 
investigate the nature of P n t, we compare the deprojected 
Pth profiles of several CF sectors with the azimuthally- 
averaged deprojected profiles of their host clusters. This 
can be done for A2204, A496, 2A0335, and Virgo, as 
shown in Fi gure El The radial means of the latter three 
clusters [32ll34j were published separately from the CF 
analyses; we divide their P t h by a fudge factor ~ 1.7 of 
an unknown origin in order to match the CF sectors away 
from the CFs; this normalization has no effect on the re- 
sults above. Also shown are the extrapolated P t h(r c f), 
assuming power-law pressure profiles. The pressure dis- 
continuities at the CFs are marked by rectangles. 

Figure [3] suggests that in Virgo and in the inner CF 
of 2A0335, Pth above the CF is comparable to the mean 
pressure at that radius, whereas below the CF, P t h (cir- 
cles) drops below its mean value (squares) . The pressure 
deviation is confined to the narrow region ~ [0.7, ]]r r f ■ 
This roughly corresponds to the shear flow layer [fj, , 
suggesting that P nt is associated with shear. 

However, in A2204 and possibly in A496, Pth in the 
CF sector deviates predominately above the azimuthal 




FIG. 3: Deprojected Pth profiles in CF sectors (circles, empty 
red or filled blue) and in their host clusters on average (green 
squares), multiplied by r to highlight the CFs. Connecting 
lines are guides to the eye. CFs are marked by rectangles 
(dashed for the new CFs, in A2204 and 2A0335). For clarity, 
Pth is multiplied by a factor of 2 in A496, and by 0.5 in 2A0335 
(see labels). Data: The mean profiles of A496 and 2A0335 are 
from [32| > Vir go f rom [34| , and A2204 from averaging the four 
quadrants in_p|. The CF profiles for A2204 are from 
A496 from 25], 2A0335 from [H and Virgo from [3J. 



average just above the CF, and returns to the mean value 
slightly below the CF. The P t h deviations are roughly in 
the region [0.9, 1.5]r c /, and are persistent although not 
highly significant. At face value, they suggest that while 
P nt just below the CF is enhanced, P n t just above these 
CFs is lower than typical of that radius. In such clusters, 
the nonthermal (and probably shear) layer may extend 
far beneath the CF, dominating much of the core. 

However, a more detailed analysis is needed to reliably 
map P n t{r). The azimuthal averages used here are sensi- 
tive to normalization problems (although not in A2204), 
non-sphericity, and deprojection errors. In addition, the 
radial mean is biased by the CFs themselves, if the un- 
derlying high density layers are spatially extended as the 
data suggest. Here, the bright layer below a CF domi- 
nates the emission at that radius, and its enhanced P n t 
may effectively lower the azimuthally averaged Pth- 

The nonthermal pressure is mostly magnetic. — The 
CF phenomenology strongly constrains the nature of the 
nonthermal pressure, as the effect is (i) robust across dif- 
ferent clusters; (ii) seen at various radii, both inside and 
outside the peak temperature radius; (Hi) is mostly lo- 
calized at CFs; and (iv) does not spread above the very 
narrow CF transitions. Shear-generated magnetic fields 
in the shear layers found beneath all CFs ll| naturally 
fulfill all these conditions. Three alternative explana- 
tions are considered in the discussion below: high-energy 
particles, small-scale turbulence (microturbulence) , and 
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TABLE I: The most significantly nonzero magnetic fields be- 
low the CF, assuming t\b — rj nt , sorted by significance. 
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a multiphase plasma. None of these are found to provide 
a viable explanation for the observations. 

We conclude that the most natural, and the only self- 
consistent, interpretation of the nonthermal pressure is 
that it is predominantly if not entirely, magnetic, in ac- 
cordance with the CF-based motivation outlined above. 
Indeed, a predominantly magnetic rut = 1 — £ — 0.2 
can explain the stability of CFs [ll, 37 1, the sharpness of 
the discontinuity, and the CF-minihalo connection. Our 
most significantly nonzero magnetic field measurements, 
assuming t/b = r)nt, are given in Table H] 

The discontinuity in the magnetic field must be in the 
tangential component, as V • B = 0, in accordance with 
the shear magnetization model Such strong fields 

are sufficiently strong to suppress the heat-flux buoyancy 
instability 38], which would otherwise tend to radially 



orient the magnetic field in the inner core. 

Discussion. — By analyzing all 18 deprojected core CF 
profiles from the literature, including three newly identi- 
fied CFs (Figured]), we discover small but significant (see 
Eq. (|TJ|) drops in P t h below the CFs. Such a discontinu- 
ity is seen in all non- merger CFs (Figure [lj , while con- 
trol samples (Figure [5]) show no such trend; systematic 
effects are ruled out. The most natural interpretation is 
an enhanced nonthermal pressure, most likely magnetic, 
confined beneath the CF. Here, the nonthermal fraction 
reaches rj nt ~ 20% of equipartition; it declines inward but 
is still noticeable ~ 1/3 and sometimes more (see Figure 
[3]) of the distance towards the center of the cluster. 

Such strong magnetization suffices to explain the CF- 
minihalo connection, stabilize the CF, and suppress 
transport across it, as needed to keep the discontinuity 
sharp. The strong magnetization, its sharp rise below 
the CF, and its coincidence with the shear layer, suggest 
saturated magnetic amplification by shear. Note that 
magnetic saturation at similar levels, i]b ~ 0.1, is some- 
times inferred from observations of radio-bright regions 
in the ICM, in radio halos and relics [3{|. The typical 
amplitudes of the magnetic field, demonstrated in Table 
HI broadly agree with Faraday rotation measurements [l9[ 
and minihalo estimates fl7j ]. 

Although our results require a predominantly magnetic 
P n t below the CF, < 1/2 of the localized nonthermal 
pressure could have a different nature. Consider first 
a high-energy particle contribution to P n t- Its pressure 



fraction, estimated as r\ ~ 10 
cool-core 17 1 and merger [3^ 
cally negligible. Shear acceleration [4l| is too slow to be 
efficient in the ICM. AGN feedback is strongly dimin- 
ished at r > few 10 kpc CFs. In conclusion, a diffuse 
high-energy component is negligible at the CFs. AGN 
bubbles could significantly contribute to P n t] here, the 
associated magnetic component must be large. 

Next, consider microturbulence. No direct measure- 
ment of turbulent pressure is presently available. Upper 
limits on the projected velocity dispersion, at the level of 
j] n t < 5-20%, were imposed in several clusters [35j, |42 - 
|44|. However, the apparent CF stability and narrowness 
imply that microturbulence cannot dominate the local 
enhancement of P nt . In particular, this strongly restricts 
the allowed level of shear-induced microturbulence. Here, 
Pb at least as strong as the turbulent pressure is required 
in order to prevent the turbulent flow from broadening 
and disrupting the CF. However, such strong, t\b > 0.1 
magnetization would stabilize the KHI, thus preventing 
the growth of turbulence in the first place. 

Deviations from a single phase ideal gas can produce 
erroneous pressure profiles, rendering our fit inappropri- 
ate (e.g., [25j). Such effects are unlikely to significantly 
distort our results, as this would require the same type of 
error over a non-negligible range of temperatures, a wide 
range of densities, and in the different environments in- 
side and outside the temperature peak. Moreover, recent 
work shows that (i) the (projected) temperature jump 
across CFs is larger (by ~ 10%) in a multiphase model 
[III; and (ii) X-ray observations of CCs are consistent 
at any given position with a single plasma phase [46j |. 
Finally, note that different multiphase plasmas on each 
side of the CF would still require strong magnetization, 
to stabilize the CF and isolate the plasma. 
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